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Adiabatic asymmetric scattering of atoms in the field of a standing wave 
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A model of the asymmetric coherent scattering process (caused by initial atomic wave-packet 
splitting in the momentum space) taking place at the large detuning and adiabatic course of 
interaction for an effective two-state system interacting with a standing wave of laser radiation 
is discussed. We show that the same form of initial wave-packet splitting may lead to different, 
in general, diffraction patterns for opposite, adiabatic and resonant, regimes of the standing-wave 
scattering. We show that the scattering of the Gaussian wave packet in the adiabatic case presents 
refraction (a limiting form of the asymmetric scattering) in contrast to the bi-refringence (the 
limiting case of the high-order narrowed scattering) occurring in the resonant scattering. 
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A. Introduction 

The observation of asymmetric diffraction in experi¬ 
ments involving the scattering of sodium atoms by a field 
of two short counterpropagating pulses of laser radiation 
m. m has stimulated several developments 12]-[I2] in¬ 
tended to explore these peculiarities of the Kapitza-Dirac 
diffraction |13j.|14j for applications in atom interferome¬ 
try [H], [H] and atom lithography m, ns using atom 
optics techniques [19], [20]. These efforts have led to 
advanced representations on the scattering of atoms by 
standing waves extending the diversity of the scenarios 
of interference occurring during the interaction of atoms 
with the field of optical lattices and, in general, mechan¬ 
ical action of light on the matter waves mm- 

The asymmetric scattering model employs secondary 
quantum-mechanical interference during interaction with 
the radiation field to achieve different intended target 
states. This interference is due to superposition initial 
states. It has been shown that the preparation of par¬ 
ticles prior to interaction in specific (in general, opto- 
mechanically mixed) states is able to dramatically alter 
the interaction pattern 0-11]. A basic example of such 
a change is the strong asymmetry in the scattering pat¬ 
tern in the case when the atomic wave packet is initially 
split into two momentum peaks differing by an odd num¬ 
ber of photon momenta [4] . Even more advanced are the 
various elaborate initial superposition states 0-0 that 
may result in a large amplitude coherent accumulation 
of the momentum on the internal energy levels caused 
by: single photon exchange [5] , narrowing of the interfer¬ 
ence fringes of the diffraction pattern [7], standing-wave 
refraction of atoms with initial Gaussian distribution of 
amplitudes by momenta etc. These effects suggest 
more flexibility in the control of atomic motion and hence 
can be useful in atom optics, in particular, in atom inter¬ 


ferometric and atom nanolithographic applications (see, 

e.g., ES]). 

The peculiarities of asymmetric scattering are ex¬ 
pressed further when dealing with the close neighborhood 
of exact resonance or when the fast switching on/off of a 
laser pulse is involved in the process. This is because, in 
these cases, stronger excitation of the system is achieved. 
Besides, the sudden inclusion of the interaction, an essen¬ 
tially non-adiabatic process, suggests more flexibility in 
choosing different preparation states. However, it is un¬ 
derstood that many of the explored effects can also be 
observed in the adiabatic regime, i.e. at large detunings 
of the wave frequency and the slow course of the inter¬ 
action. Since the adiabatic interaction schemes as a rule 
suggest more robust technologies, complementary discus¬ 
sions of the adiabatic model of the asymmetric scattering 
are in demand to clarify the potential for controlling the 
diffraction picture of atomic wave-packets in this regime. 

In the current paper, we first present the simplest 
scheme for the coherent diffraction of atoms by a standing 
wave occuring strongly asymmetrically in the adiabatic 
regime - the case where the initial atomic wave packet 
involves only two translation states whose momenta dif¬ 
fer by two photon momenta. We then reveal that the 
same form of initial coherent superposition state of the 
atom may cause qualitatively different effects for adia¬ 
batic and resonant standing-wave scattering (though this 
is not necessarily the case for all initial wave-packets). 
Finally, we demonstrate that the evolution of the Gaus¬ 
sian momentum distribution of amplitudes at adiabatic 
standing-wave scattering presents a strongly asymmetric, 
with respect to the initial momentum direction, diffrac¬ 
tion (with minor deformation of the wave packet form) 
while at resonance the same wave packet undergoes sym¬ 
metric scattering (analogous to the high order narrowed 
diffraction discussed in 0)- 
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It should be noted that the scattering model in the 
form presented below covers only some of the important 
features of the asymmetric diffraction. For instance, it 
does not describe the oscillatory dependence of the scat¬ 
tering amplitude on the resonance detuning [2]. How¬ 
ever, the inspection of the known developments for the 
resonant case shows that it is possible to include these 
peculiarities by means of modification of the intermediate 
preparation states and taking into account the first-order 
non-adiabatic corrections. 


B. Adiabatic asymmetric scattering 


The dynamics of an effective two-state system in the 
field of a standing wave E = 2i?o/(t) cos{kz) cos(wt) with 
a slowly varying envelope f{t) at small interaction times 
(in the absence of spontaneous emission) is described, in 
the rotating wave and Raman-Nath approximations, by 
the time dependent Schrodinger equations for the prob¬ 
ability amplitudes of the states 01^2 


= 211^1 {t)cos{kz)e *^^* 02 , 
= ‘2Uof{t) cos{kz)e~^"^^ai, 


( 1 ) 

( 2 ) 


where A = 0 J 21 — wg is the resonance detuning, Uq = 
—dEo/(2h) is the peak Rabi frequency of the travelling 
wave, d is the dipole moment of the transition under 
consideration. 

In the adiabatic regime of large detuning, slow inclu¬ 
sion and variation of the interaction, |A| t ^ 1, this sys¬ 
tem is reduced, via adiabatic elimination of the excited 
state (see, e.g., m\ ,[28]), to a simple first order equation 


i 


dai 

dt 


4|I/o|^ 

A 


f^{t) cos'^{kz) ai, 


(3) 


we suppose that before the preparation in the state ([^ 
the atom had an exactly defined momentum pq, that is 
/ p{z)e-'^^^/^d{z) = 6{p-po)): 


ai = e 


■• 2|Uq| 


^ ^(signA • i) 


..w 2 i + (-ir 


^^(signA-i) S{p - po - nhl^] 

where J is the Bessel function. The corresponding prob¬ 
ability of absorbing n photons is written as: 


W^{t) = 


1 • A ■\-‘m r /2|f7o|^ 

-^- ^(signA • l) a2mJnl2-m I 

m \ 

( 8 ) 

If the initial wave packet is not split in the momentum 
space, i.e. if a 2 m = 0 at m 7 ^ 0, then the solution Q 
becomes the well-known expression mi,mi: 


VF„(t) = 


l + (-l )%2 2 |t/o|"r 


4/2(w): W = 


|A| 


(9) 


that describes symmetric, with respect to the initial 
atomic momentum po, diffraction pattern. 

However, as it is readily seen from Eq. (|^, the situ¬ 
ation is qualitatively changed at the initial conditions of 
splitting. It is seen that a secondary interference takes 
place; during diffraction by the standing wave the diffrac¬ 
tion peak sets {Jn/ 2 -m}! originated from the correspond¬ 
ing peaks of the initial wave packet, overlap. This inter¬ 
ference significantly changes the scattering pattern. In¬ 
deed, consider, for instance, the simplest case when the 
initial wave packet is split into only two peaks: 

ai( 0 ) = (ao + 0 : 26 *^'''') 

02 ( 0 ) = 0 . ( 10 ) 


the solution of which is straightforward: 

. 2 \ Uo \^ . 2 \ Ur )\^ r ,, 

ai(t) = ( 4 ) 

where r is the integral of the square of the field envelope: 


I 

T = J 


(5) 


For the probability of the n-th diffraction order we then 
get a strongly asymmetric scattering: 


Wnit) = 


! + (-!)” 


(X2 


o^oJn/ 2 {u) + j signA 


( 11 ) 

Indeed, taking into account the equality J_„ = (—1)" J„, 
we have: 


0 

To present the simplest model for the asymmetric scat¬ 
tering in the adiabatic regime, we consider, following [7], 
the scattering caused by the initial conditions of the form 

+ 00 

ai( 0 )= ^ a 2 „e* 2 -fcXz), 

m——oo 

02 ( 0 ) = 0 , ( 6 ) 

which, evidently, can be created using adiabatic pro¬ 
cesses. The solution of the coherent diffraction prob¬ 
lem in the momentum representation is (for simplicity. 


Wnit) = 


l + (-l)^ 


dn/2 + I« 2 |^T 2 / 2 - 1 - 

Jn/2Jn/2-l, n > 0, 

l“ol •^^n/2 + l“2|^ -l^„/2+i + 

d-n/2J-n/2+l, n < 0. 

^ ( 12 ) 

We note that the scattering probability (11) has the 
same structure as the corresponding non-adiabatic prob¬ 
ability for the exact resonance case with preliminary ex¬ 


citation of the atom by a travelling wave |4]. Hence, the 
peculiarities of the diffraction process in the adiabatic 
and non-adiabatic regimes are qualitatively the same. 
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The scattering pattern asymmetry is defined as 
+ 00 

AW{t) = Y,iW+n-W_^) 

n—1 

1 ^ |2 / 72 I r 2 \ ^ Im(Q^oQi 2 ) 7 ^ r T \ 

— \o-2\ {^0 -fI) - signA — ~ 

u 

Co = JiJoiu) + JUu))du, Co |„^oo « 0.63^13) 
0 

Consequently, the maximum possible asymmetry is 
achieved at 


Im(aoQ!2) = ±1/2 


(14) 


(|ao| = \cx. 2 \ = l/-\/2); and it is the same as in the non- 
adiabatic case: at u —>■ 00 , more than 80% of atoms 
deflects to a definite direction. However, the preferable 
direction of the deflection is determined, as it is seen from 
Eq. (12), not only by the sign of Im(Q;oQ: 2 )) but also by 
the sign of the detuning. The total acquired momentum 
of the atom after the interaction is 


(p) = 2nhkW2n = ‘2hk ^|q; 2|^ - Im(Q;oQf2) ^ ■ 

(15) 

The first term in the right-hand side of this equation, 
which is the momentum shift coming from the splitting 
of the initial state §, is not more than two photon mo¬ 
menta, but the second term (the result of the standing 
wave action) is not restricted. However, we note that the 
mean momentum increase is determined, in addition to 
the factor Im(Q:oa 2 ) coming from the initial conditions, 
by the parameter u = 2 \Uo\^t/\A\ which is supposed 
to be not too large within the Raman-Nath approxima¬ 
tion. The approximation is valid if the gained kinetic en¬ 
ergy of the atom Ckin = {nhk)'^/(2m) is small compared 
with the interaction energy e = hUo for all diffraction 
orders n. It follows from the properties of the involved 
Bessel functions that the maximum populated diffrac¬ 
tion order is approximately Umax ~ 2u. Hence, should 
be Au^LOrec ^ min/I/o, 1/t}, where w^ec = kk^/(2m) is 
the recoil frequency and t is the interaction time. The 
diffraction pattern caused by the initial splitting ( [Io| ) at 
the maximum possible asymmetry is presented in Figjl] 
For the chosen parameters Uq and A this pattern is con¬ 
sistent with the Raman-Nath approximation, e.g. for 
sodium atoms with m = 23 amu and optical field with 
A « 0.5 /rm, for interaction times r < 10“^ s. 

Thus, we have seen that the adiabatic interaction with 
the standing wave can also occur asymmetrically pro¬ 
vided the initial wave packet is split in a special man¬ 
ner. It should be noted that, interestingly, the wave- 
packet ( |10[ ) has the same form as the one considered in 
m (Eq- (21)) when discussing the narrowing of the inter¬ 
ference fringes of the diffraction pattern at the exact res¬ 
onance. Thus, the same form of initial wave-packet split¬ 
ting may, in general, lead to different diffraction patterns 
for adiabatic and resonant standing-wave scattering. 



FIG. 1. Diffraction patterns at adiabatic interaction of atoms 
with a standing wave: (a) usua l deflection - Eq. (b) 

asymmetric deflection - Eq. (Ill, ao = 02 = 

At = —500, Ut = 50. 


Nevertheless, the wave packet ( |10[ ) itself can not be 
used to achieve both resonant narrowed (at |Ar| «C 1) 
and adiabatic asymmetric (at |Ar| 1) scatterings. 
The narrowing condition a 2 = —ao [Z] and the condi¬ 
tion for maximum asymmetry (14) are not compatible 
since the phase conditions that should be imposed on 
ao and 02 contradict. However, this incompatibility is 
not necessarily the case for all the possible wave-packets. 
We demonstrate this below by examining the behavior 
of the Gaussian wave packets. We will see that the adia¬ 
batic standing-wave scattering of a Gaussian wave packet 
presents refraction (a limiting form of asymmetric scat¬ 
tering, see also m) in contrast to the bi-ref ring ence (the 
limiting case of high-order narrowed scattering, H) oc¬ 
curring with the same wave-packet at the resonant scat¬ 
tering. 


Regarding the preparation of atoms in the states §, 
first we note that such distributions can be achieved, in 
general, in effective two-state systems when the terms 
a 2 raC‘^™‘^^^{z) correspond to the different levels that 
compose the effective ground state. Such a situation is 
the case, e.g. in the two-level systems with magnetic sub- 
levels (see an example of such a preparation of atoms by a 
single elliptically polarized travelling wave pulse in m)- 
Alternatively, this distribution can be viewed as one cor¬ 
responding to the same internal state of an atom. For 
instance, a combination of adiabatic rapid passage and 
multiphoton Bragg diffraction can be used to efficiently 
transfer many photon recoils of momentum hk to cold- 
atoms, thus creating two- or multi-peak distributions of 
the needed structure [3T]. Creation of a Gaussian mo¬ 
mentum profile starting from the atoms prepared in a 
very narrow distribution around p — 0, then adiabati- 
cally ramping on a one-dimensional optical lattice and 
further suddenly spatially shifting the lattice by 1/4 of 
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FIG. 2. Representations of the same Gaussian initial wave 
packet (M = 10) for adiabatic (a 2 m) and resonant(sm) scat¬ 
tering regimes. 


the lattice period, is reported in [32]. 

Further, as a systematic method to create the desired 
distributions of the populations in the multi-level systems 
one may apply different STIRAP schemes |29| involving 
combinations of traveling and standing waves. An exam¬ 
ple of such a process leading to the creation of discrete 
wave-packets of exponential distribution by momenta is 
suggested in |7]. Momentum state preparation of a two- 
level atom using two-stage Kapitza-Dirac diffraction of 
an initially single-momentum atomic beam was recently 
discussed in |33| . Gaussian wave packets can be prepared 
by exciting atoms to Rydberg states [34] (for preparation 
of molecular wave packets by femtosecond pulse technol¬ 
ogy see, for instance, [55]). 


C. Diffraction in the case of initial Gaussian 
distribution by momenta 


For a discrete Gaussian wave-packet involving only 
even orders (see Figl§: 

gi(a-|-(signA)'n-/2)m 2 

(^m)V4 (16) 

where M is the distribution half-width, the adiabatic 
scattering probability ^ is rewritten as (i^ = n/2 — m) 

2 


W„. = 


_ i+(-iy 

2 


■>/2 


+ 00 




E e 

V — — 00 


— lOLV—ynl — 






(17) 


where u = 2 \Uq\ r/|A|. The behavior of the functions 
In{u) was studied in m- The approach employs the 
following exact linear differential-difference equation: 


din _ n u , 


2M 


^ T 

Jrj.— 1 


(18) 


which is readily derived by differentiating and using 
the identity 2vJ^{u) = -\- Jv-i)- 

Examining the solution of this equation, we note that 
at the beginning of the scattering, u = 0, the solution is, 
as expected, the Gaussian 

n2 

In = coe~^ (19) 

with Co being the pre-factor of the exponent in Eq. ([T^ . 
Eurthermore, dividing the equation by n/M and passing 
to the variable I = n^/(2M) (<J4> dl = {n/M)dn) we see 
that the term in the brackets in Eq. ( [T^ becomes pro¬ 
portional to u/(2n), hence, it can be neglected for large 
diffraction orders such that \n\ 3> u. Thus, the dynamics 
n of the wave packed given by /„/ 2 (m) is effectively local¬ 
ized within the interval —2u < n < 2n. We will see 
that during that time the wave-packet moves as a whole 
within this interval and slightly changes its form. 

To discuss the diffraction details, we expand In±i into 
the Taylor series at the point n and keep the first three 
terms. The resultant equation reads 


din 

dn 



u 

M 


(cos a) 


In + 


I d^In 
2 dn^ 


+ l(sinQ:) 


din 

dn 


( 20 ) 

If cos a = 0, the equation is reduced to a first-order 
one, the solution of which is again given by a Gaussian: 


In = Coe , ( 21 ) 


where now cq = [ttM (l-|-The distri¬ 
bution half-width is hence, the wave-packet 

steadily broadens during the time. 

If cos a ^ 0, equation (201 is reduced to the Airy 


equation m- Accordingly, the solution of the diffrac¬ 
tion problem finite at n —>■ ±oo in this case is written in 
terms of the Airy function of the first kind |5S|. How¬ 
ever, before discussing this solution, it is helpful to take 
a look at the solution derived if one only neglects the 
second-derivative term in Eq. (20): 


(,.-(coso)u )2 

In = Coe 2 (ai+.(=.„o)„ 


( 22 ) 


where now cq = [ttM (l-|-(sina)^u^/M^)] . This 

solution indicates that in the case of non-zero cos a the 
wave-packet in general broadens and moves as a whole in 
the momentum space, the broadening being defined by 
sin a and the displacement being proportional to cos a. 

Consider now the exact solution of Eq. (201: 


4 = Coe'^"Az(A + h^), 


(23) 


where the parameters A, h are written as 
^ ^ n - (cos q)u ^ ^ M -i{s\na)u 

2“i/3 (|cosq:| u)^^^ 21/3 (|cosq;| 

and the constant cq is defined from the normalization 
condition. This is a localized wave-packet the proper¬ 
ties of which are very controlled by the imaginary part 
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of the argument z = iV + of the Airy function: 
Im( 2 :) = Im(/i^) ^ sin(Q;). Indeed, consider the case 
sin a = 0 when the argument is real. Then, arg 2 = 0 
for n > u and arg z = —tt for n < u. The asymp¬ 
totes of the Airy function for large real argument z = x 
are known to be Ai{x) ~ jx^/‘^ if a: > 1 and 

Ai{x) ~ sin(a: -I- 7 r/ 4 )/a:^/"^ if a: < —1 |3n]- Hence, the 
localization at the packet side for which iV < 0 is due 
to the factor (since Re(/i) is positive) and the wave 
packet is localized at its other side due to the Airy func¬ 
tion asymptote. 

If sin a ^ 0, then argz ^ 0 , 7 r for any n (because 
N is real and Im(/i) ^ 0). The real part of the ar¬ 
gument of the Airy function becomes zero near the 
point n ~ (cosa)u -I- 0(I/n). Going far away from 
this point, that is at (n — (cosa)u) —>■ ±oo, depend¬ 
ing on the sign of (sina)/u, argz either tends to zero 
or — TT. The asymptote of the Airy function for large 
\z\ is Ai{z) ^ j in the sector |argz| < tt 

and A*(z) ~ in the sector [tt — arg(z)| < e <C 1 . 

Accordingly, for the side of the distribution for which 
(sin Q:)/n is positive (Im(z) < 0 ) the localization is due 
to the factor (the Airy function adds negligibly small 
oscillations) and the localization is due to the asymptote 
of the Airy function in the opposite side of the packet. 
We conclude by noting that the Airy-function solution 
(23 1 -(24 1 presents a highly accurate approximation. The 


comparison with the exact numerical result and the mov¬ 


ing Gaussian approximation (22), which provides quali¬ 
tatively rather a good description, is shown on Figj^. 

Thus, during the time evolution, the wave packet al¬ 
ways remains localized moving in the momentum space 
if cos a ^ 0 and broadening if sin a ^ 0. If a = 7 r /2 ± 
TTk,k = 0 , 1 , 2 ... the distribution peak does not move; 
it just broadens. In contrast, if a = ± 7 rfc,fc = 0,1,2... 
the broadening is absent in the first approximation and 
the distribution displaces as a whole (see Fig[^. The 
displacement of the peak position as well as the dis¬ 
tribution broadening are determined by the parameter 
u = 2|t/o|^T/|A|. Thus, the adiabatic scattering of the 
Gaussian wave-packet (16) presents, approximately, a re¬ 


fraction to a definite angle controlled by the interaction 
time with the light field. This picture significantly differs 
from that for the resonant scattering regime. 

Indeed, consider the resonant diffraction of the same 
Gaussian wave packet ( [T6| ). The scattering probability 
in this case is written as [TT] 


w:{t) = 






(25) 


where t is the interaction time and the corresponding 
initial-state vector for a distribution of the form (§ 
is defined as 

1 (-T)"* 

Sni{rn) = - ^(26) 

(see Fig. 2). 





FIG. 3. Adiabatic diffraction pattern of the Gaussian wave 
packet with M = 10, a = tt: At = 500, Ut = 50. Dashed line 
- moving Gaussian approximation Eq.(22). Solid line - Airy 
function solution Eqs. (23),(24). 


Rewriting i^Sm as {m ^ n — v) 


g 2 Ml 




^ — 


-h (-l)”-'^e ^ 


— — — — \-i/3{n — iz) 




with Ml = 4M and 

1 


/3 = -(2a -I- 27r -I- TT sign(A)), 


(27) 

(28) 


thus splitting the sum in Eq. (25) into two parts, then 


passing to a new summation variable v = n — m and 
further using the identity {—\Y J^{u) = Ji,[—u), we get 


w: = 


1 


2i/7r 


lI„(Ur) + (-1)M„(-Ur)r 


(29) 


where = 2Ut and 
+00 


^;k) = E' 


—i^v—{n—vY/2M\ 


Ml/4 


J„(Ur). (30) 


As it is immediately seen, equation (29) describes a sym¬ 


metric two-fringe scattering pattern. Hence, the resonant 
scattering of the Gaussian wave packet (16) presents bi¬ 
refringence. This is demonstrated in Fig|4pwhere the pre¬ 
sented graph has been calculated using the exact equa¬ 
tion (25). 

We note that /(((ur) given by Eq.(30) is exactly the 


same function as Iniu) (with parameters altered as u —)■ 
Ur, M —>• Ml/A, and a —>■ /3 according to Eq. (28)) 
that we used in treating the adiabatic scattering. Thus, 
one may use the above moving Gaussian approximation 
(22) or the Airy-function solution (23)-(24) to accurately 


explore the interference of the fringes at resonant scat¬ 
tering. At a = ± 7 rfc, fc = 0, 1 , 2 ... and large enough Ur 
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D. Summary 

Thus, we have presented a model of adiabatic asym¬ 
metric scattering of coherent superposition states of ef¬ 
fective two-state atoms in the held of a standing wave. 
We have shown that the behavior of the atomic wave- 
packets at adiabatic diffraction may signihcantly differ 
from the diffraction at resonant scattering. 

We have demonstrated that the discrete Gaussian 
wave packet at adiabatic diffraction undergoes highly 
asymmetric scattering (refraction) while the evolution 
of the same wave packet at resonant scattering can 
be characterized as a high-order narrowed scatter¬ 
ing (bi-refringence). Hopefully, these peculiarities of 
the standing wave diffraction of atoms with Gaus¬ 
sian initial momentum distribution will be useful for 
atom interferometric and atom lithographic applications. 


FIG. 4. Resonant diffraction pattern of the Gaussian wave 
packet (161 with M = 10, a = n: A = 0,Ut = 50. 


when the functions In{ur) and In(—uf) practically do not 
overlap (i.e., at Ur > M), we have 

Wf « ^ {\IuiUr)\^ + |/n(-tX.)|") , (31) 

which describes a two-peak diffraction pattern, each of 
the peaks being slightly different from the initial Gaus¬ 
sian form. 
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